Field-driven transition in an Ising magnet with mixed interactions 
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Transfer-matrix methods are used, in conjunction with finite-size scaling and conformal invariance 
concepts, to generate an accurate phase diagram for a two-dimensional square-lattice Ising spin-1/2 
magnet, with couplings which are positive along one coordinate axis, and negative along the other, 
in a uniform external field. Our results indicate that the critical line starts horizontally at the zero- 
temperature end of the phase boundary, at variance with the reentrant behavior predicted in some 
earlier studies. Estimates of the thermal scaling exponent are very close to the Ising value j/t = 1 
along the critical line, except near T = where strong crossover effects prevent a reliable analysis. 
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I, INTRODUCTION 

In this paper we investigate a square-latice Ising spin- 
1/2 system with ferro- and antiferromagnetic interac- 
tions, in the presence of a uniform magnetic field. The 
Hamiltonian is given by: 



n 



(1) 

where J x , J y > 0, and <Jij = ±1. Here, all fields, cou- 
pling strengths and temperatures will be given in units 
of J x - At T = 0, H < 2J y , the ground state con- 
sists of alternating ± stripes along the x axis, while for 
larger H all spins are parallel to the field. A critical 
line T C (H) connects (T = 0, H — 2J y ) to the zero- 
field Onsager critical point, which for J y = 1 is at 
T c ° = 2/ ln(l + V2) = 2.2691853 .... Since the ground 
state does not exhibit macroscopic residual entropy, the 
transition along T C (H) is expected [l|| to belong to the 
Ising universality class for all T ^ 0, similarly to the 
closely-connected case of the standard (isotropic) anti- 
ferromagnet in a uniform field 0, H, 0, H, 0| ■ 

The problem described by Eq {!]) was treated by 
Miiller-Hartmann and Zittartz's interface method [2] in 
Ref. 0, which provides an excellent summary of earlier 
work. It was revisited in Ref. [§, using an approach which 
considers the zeros of the partition function on an ele- 
mentary lattice cycle, and their connection to the free 
energy singularity at the transition Remarkably, the 
critical line found in Ref. d is predicted to display a pos- 
itive slope close to T = 0, so the critical field reaches 
a maximum at some nonzero T before approaching zero 
at higher temperatures. Similar reentrant behavior was 
predicted upon application of a Bethe-Peierls approxima- 
tion K\. On the other hand, the linear-chain approxima- 
tion [l| gives an exponentially vanishing positive slope at 
T = 0, while the interface method also predicts an expo- 
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nentially vanishing value, only on the negative side [3]. 
Recent real-space renormalization results [lOj point to a 
finite negative slope at T — 0. 

We use transfer-matrix (TM) methods, in connection 
with finite-size scaling and conformal invariance ideas, in 
order to produce a numerically accurate phase diagram 
for this problem. The underlying hypotheses in our work 
are: (i) that the phase transition is second-order all along 
the critical line, and (ii) that it belongs to the Ising uni- 
versality class. Both assumptions are critically reviewed 
towards the end of the paper, in light of the numerical 
results obtained while assuming their validity. 

In Sec. [II] we recall the calculational methods used for 
the approximate location of the critical line; the respec- 
tive results are exhibited, as well as their extrapolation 
in the thermodynamic limit. In Sec. IIIIl we analyze the 
data generated in Sec. HU both in comparison with the 
existing literature, and in regard to their internal consis- 
tency. The universality of critical behavior is discussed, 
and concluding remarks are made. 



II. METHOD AND RESULTS 

We have kept J y — 1 in all calculations reported here. 
We set up the TM on strips of width N sites, with peri- 
odic boundary conditions across. Referring to Eq. |T|. 
three choices of orientation are available in this case, 
namely the TM can be iterated along the ferromagnetic 
(F), or x- direction; along the antiferromagnetic (AF), 
or y- direction; or along the diagonal (D) of the square 
lattice. As this is a weakly anisotropic system [ill. fl2]|. 
one would expect estimates of, e.g., critical exponents 
and locations of critical points, to converge to the same 
orientation-independent limit for N 3> 1, while finite-size 
corrections should differ in each case. In order to obey 
the ground-state symmetry, only even values of N are 
allowed for F and D, while no such restriction applies 
to AF. We generally used 4 < N < 20; for F we went 
up to TV = 22. This range of N enabled the authors 
of Ref. 5 to locate the critical line of the isotropic Ising 
antiferromagnet to very high accuracy. 
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Figure 1: (Color online) Low-temperature approximate crit- 
ical boundaries given by solutions of Eq. ([2]). Triangles, 
squares, and hexagons denote respectively TM along F, AF, 
and D direction. Empty symbols: N = 4; full symbols: 
N = 22 (F), 20 (D and AF). 

A. Keeping 77 = 1/4 

Following earlier work on similar problems 0, 0, [T3 |. 
our finite- TV estimates for the critical line are found by 
requiring that the amplitude-exponent relation of confor- 
mal invariance on strips [3] be satisfied, with the Ising 
decay-of-correlations exponent r\ — 1/4: 

4Nk n (T,H)=tt , (2) 

where kn{T,H) = In | A1/A2 | is the inverse correlation 
length on a strip of width N sites, with Ai, A2 being the 
two largest eigenvalues (in absolute value) of the TM. 

For H <C 1, the solutions to Eq. leave the T axis 
vertically, and are very close to each other for all three 
orientations of the TM. We shall briefly return to this, 
towards the end of the paper. For low temperatures T < 
1, the main region of interest here, substantial differences 
arise. These are illustrated in Fig. [TJ One sees that, 
with increasing N, all three families of curves get closer 
together, though AF does so at a slower rate. Also, in all 
cases their limiting shape for T < 0.4 is, with very good 
accuracy, a straight line going through the exact zero- 
temperature fixed point at H = 2. Reentrant behavior is 
clearly visible for the F curves and, to a lesser extent, for 
the D ones, though in both cases the peaks turn flatter as 
N increases. In order to check whether the reentrancies 
vanish in the N — ► 00 limit, we plotted the sequences of 
slopes Sn of the T — > straight-line sections referred to 
above, against N~ x . 

Results are shown in Fig. [2j For both F and D the 
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Figure 2: (Color online) Slopes Sn of the low-temperature 
straight-line sections of approximate critical boundaries given 
by solutions of Eq. |[2|), against 1/N. Triangles, squares, and 
hexagons denote respectively TM along F, AF, and D direc- 
tion. For AF, data are scaled by a factor of 5. Lines are fits 
to data: linear for F and D (all N), and parabolic (linear plus 
quadratic) for AF (N > 10 only). 



trend is, to an excellent approximation, Sn — > linearly 
with TV -1 . Indeed, the respective extrapolated values 
are (4 ± 3) x 1CT 7 (F) and (-15 ± 4) x 1CT 5 (D), which 
in practice equate to zero within the present context. 
For AF, a large amount of curvature is present; from a 
parabolic fit (S N = S ao + aN- 1 + b N~ 2 ) of N > 10 data, 
one gets Soc = —0.029 ±0.004. This would correspond to 
1 - H c (T)/H c (0) ~ 0.01 at T = 0.5 which, though small, 
is significant. 

In summary, the analysis of limiting slopes at low T, 
N — ► 00, indicates that the reentrant behavior observed 
for F and D data is a finite-size effect, which tends to 
vanish in the thermodynamic limit. A small discrepancy 
still remains, between F and D data which are consistent 
with a horizontal critical line at T — > 0, and AF results 
which point to a slightly negative slope in that limit. 



B. Phenomenological Renormalization 

Further insight can be gained by relaxing the assump- 
tion made in Eq ^ , that the phase transition belongs to 
the Ising universality class, and demanding only that it 
remain of second order. From finite-size scaling, one gets 
the basic equation of the phenomenological renormaliza- 
tion group (PRG) for the critical line: 

Nk n (T,H) = N' Kn ,(T,H), (3) 
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Figure 3: (Color online) Low-temperature approximate criti- 
cal boundaries given by solutions of Eq. ([3]). Triangles: TM 
along F direction, TV = 14; hexagons: D, TV = 14; empty 
squares, crosses, and full squares all for AF, respectively 
TV = 10, 12, and 14. 



Figure 4: (Color online) Slopes Sn of the low-temperature 
straight-line sections of approximate critical boundaries given 
by the solution of Eq. for TM in the AF direction, against 
1/7V. Data points are for TV = 5 — 14. The line is a parabolic 
(linear plus quadratic) fit to data. 



where the strip widths TV and TV' are to be taken as close 
as possible for improved convergence of results against 
increasing TV. This means TV' = TV — 2 for F, D, and 
TV - 1 for AF. 

We found that PRG results converge very rapidly for 
both F and D orientations of the TM, without any sign 
of the reentrances shown by the solutions of Eq. J2]) . For 
the F case, the largest discrepancies between TV = 6 and 
8 amount to 0.8% close to H — 0, and are slowly reduced 
upon increasing H; around T ~ 0.7, H ~ 1.97, the two 
curves differ by 0.04%. At T = 0.4, both coincide to 
within two parts in 10 5 , at less than 0.1% of H = 2, 
and then home in towards (T, H) = (0, 2) on a straight 
line. The discrepancy between TV = 12 and 14 is never 
more than one part in 10 4 for T < 1. The picture is 
quantitatively similar for the D orientation. 

On the other hand, for PRG with the TM along the AF 
direction, one gets relatively large negative slopes (but 
approaching zero with increasing TV) as T — > 0. Further- 
more, at intermediate temperatures 0.4 < T < 0.8 the 
curves show inflection points which make extrapolation 
of such sections to TV — > oo prone to instabilities. 

The above-mentioned features are illustrated in Fig. [3j 
We investigated the behavior of the limiting slopes for the 
AF family against 1/TV, with results displayed in Fig. [H 
The line shown in the Figure is a parabolic fit to data 
(Sn — S'oo+a TV -1 + b TV -2 ), from which one gets \Soo\ < 
10 -3 , i.e., essentially zero. 



C. Extrapolations 

In this Subsection, we deal directly with extrapolations 
of finite-size data for all (T, H ) . This is in contrast to the 
analysis of slopes, which applies only at low T where the 
finite- TV data actually fall on straight lines. 

Extensive investigation of the related problem of 
isotropic Ising antiferromagnets in a field [Sj shows that, 
in that case, the main irrelevant exponent is j/j r = —2, 
i.e., the leading corrections to scaling are expected to de- 
pend on TV -2 . Thus, it is plausible to assume that such 
corrections also play a dominant role here. However, the 
analysis of slopes above suggests that corrections in TV" 1 
are present as well, at least in the low-temperature re- 
gion. 

In what follows, we examine both scenarios, i.e., yi r 
1 and —2. The locations of points on the approximate 
critical line, for strip width TV, are denoted by (T N , H N ). 
On account of the overall shape of the phase diagram, we 
extrapolate against negative powers of TV in two different 
ways: (i) at constant X = H for high T > 1.7, and (ii) 
at constant X — T for T < 1.7. For large TV, we fit Y/v 
[= T N {H) in (i), or H* N (T) in (ii)] to a form 

Y N (X) =Y cxt (X) + a yir (X)Ny- +b Vii (X)N 2 y- . (4) 

We have chosen to calculate Y cxt (X), a Vil (X), and 
b yil (X) from the three largest values of TV available for 
each strip orientation, thus error bars (other than those 
associated to the positions Y N (X) themselves) are not 
available from this procedure. 
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Figure 5: (Color online) Low-temperature extrapolated crit- 
ical boundaries, obtained via Eq. Q, with yi v = —2 (empty 
symbols), and j/i r = — 1 (full symbols). TM is in the F direc- 
tion in (a), and along the diagonal (D) in (b). 



Figure 6: (Color online) Low-temperature behavior of fit- 
ting coefficients a_i (empty symbols), b_i (full symbols) of 
Eq. (H}. Triangles: TM along F direction. Hexagons: TM 
along diagonal (D). Inset: high-temperature behavior. Same 
axis labels and symbol captions as in main figure (see text). 



We managed to produce well-behaved extrapolated 
critical lines from the solutions to Eq. ([2]) , only for F and 
D orientations of the TM. For AF, though extrapolations 
are generally smooth for T > 0.6, they display instabil- 
ities for lower T, i.e. the section of the phase diagram 
which is most relevant in the search for reentrant behav- 
ior. As regards PRG curves, due to the fast convergence 
of F and D data we found that N = 14 data can already 
be taken as very close to the N — > oo limit, to within 
an estimated one part in 10 4 . As mentioned above, the 
extrapolated PRG curves for AF display instabilities, ex- 
cept for very low T (for which the relevant information 
is summarized in the slope analysis illustrated in Fig. [4j) . 

Fig. [5] shows the low-temperature regions of extrapo- 
lated critical boundaries, for F and D, assuming y lt = — 1 
or —2 in Eq. - For both curves corresponding to 
y- n = — 1, there is a rather flat section: for F, points with 
T < 0.36 remain within one part in 10 5 from H = 2, 
while for D the T < 0.35 region is within one part in 
10 4 from that limit. On the other hand, the y lr = —2 
extrapolations exhibit tiny reentrances, with maxima re- 
spectively at H = 2.0024 (F) and H = 2.0012 (D). 

In order to get further insight into the competing sce- 
narios under investigation, we examine the behavior of 
the coefficients a_i and 6_i of Eq. ([4]) along the critical 
line. Indeed, this amounts to an unbiased test of whether 
the dominant N- dependence of finite-size data is on A -1 
(via a_i) or N~ 2 (via In Figure[6]one sees that, at 

low T, the A -2 terms tend to vanish for both F and D 
cases. Furthermore, the inset of the Figure shows that in 
the low-field limit, it is the N~ 2 corrections that become 



dominant, and those in A -1 become negligible as H — > 0, 
as is well known for the zero-field Ising model 0, EH] • The 
latter result gives further credence to the procedure just 
described, and thus to the conclusions reached regarding 
the low-T regime. 

Finally, we found no clear evidence that an exponent 
2/vac = —4/3, associated to vacancy excitations, might 
be present (3|]. This is because (i) at the F and D data 
clearly depend predominantly on A" 1 at low T, as al- 
ready shown in the slope analysis; and (ii) the slope data 
for AF do not exhibit any significant improvement in 
quality of fit when plotted against A~ 4 / 3 (compared to 
the parabolic fits in A -1 of Figs. [2] and |4]). 



III. DISCUSSION AND CONCLUSIONS 

We begin by recalling that Eq. J3]) depends on the exis- 
tence of an underlying diverging correlation length, thus 
its non-trivial fixed points correspond to a second-order 
transition. We only failed to find such points for low 
temperatures, generally T < 0.2. This is due to the ex- 
tremely slow convergence of numerical results, connected 
to the very large ratio between positive and negative ex- 
ponentials which are the TM states' Boltzmann weights 
in that region. If a tricritical point were present, separat- 
ing first- and second-order sections of the critical curve, 
one would expect spurious effects such as the "hooking" 
found for the three-dimensional version of the current 
problem On the contrary, as long as we can find 
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low- T solutions to Eq. ((3j) they behave in the expected 
manner, i.e. homing in towards the exact T = fixed 
point at H — 2. Furthermore, the extrapolated r\ = 1/4 
curves agree very well with the solutions of Eq. ([3]) down 
to T = 0.2, and still extend somewhat further down to 
T sa 0.1 Since, by conformal invariance, rj — 1/4 corre- 
sponds to an Ising transition, our results indicate that 
this is the character of the critical line, at least down 
to T = 0.1. Therefore, if a tricritical point is present, 
it must be located at T < 0.1, H » 2. We thus con- 
clude that the transition is indeed second-order and in 
the Ising universality class along the whole of the critical 
curve (except for the latter possibility, which we are not 
able to probe directly). 

We now recall that the reentrant behavior predicted in 
Ref.H is sizable: in the units used in the present work, it 
translates into the critical line leaving (T, H) — (0, 2) 
with a slope S = (1/2) In 2 = 0.3466... [see their 
Eq. (31)]. This is in contrast with the results of Sub- 
sec. [TIC] where we find at most S = 7x 10~ 3 (with the 
TM in the F direction, using y„ = —2). 

In the comparable problem of isotropic antiferromag- 
nets, though the critical line H C (T) given in Ref. @ does 
not exhibit reentrances, it is always above that found in 
Refs. |4j, |5j (except at the T — and H — ends, where 
both lines coincide). The maximum discrepancy, of or- 
der 4%, is in the central region, 0.7 < T < 1.5, tailing 
off towards both ends. It thus appears that the methods 
employed in Refs. @, d generally tend to overestimate the 
extent of the ordered region in parameter space. 

Though the present problem is weakly anisotropic in a 
broad sense, the distinct nature of spin couplings along 
each coordinate axis is responsible for the introduction 
of subtle biases, when one iterates the TM along ei- 
ther of those very same axes. Indeed, the size of the 
ordered region predicted by low-T TM results system- 
atically decreases as one changes orientation from F to 
D, and finally to AF. The explanation is that, for the 
high fields near the critical curve, the ferromagnetic cor- 
relations picked out when the TM goes along F tend to 
be emphasized; when the TM goes along AF, the corre- 
sponding antiferromagnetic correlations are inhibited by 
the field. The small negative extrapolated slope of the 
AF curves, shown in Fig. most likely reflects the lat- 
ter effect. For the TM going along D, an evenly-balanced 
mixture of both kinds of correlations is collected upon its 
iteration. Therefore we believe that, of the three setups 
implemented here, the latter is the one likely to produce 
the most reliable results. 

Note also that, for PRG one is always comparing corre- 
lation lengths evaluated along the same lattice direction 
in Eq. ([3]), so the aforementioned biases tend to cancel 
out, if present. Thus, (i) the PRG curves for F and D 
do not show reentrancies, even for finite N; and (ii) the 
slope of PRG curves for AF approaches zero as N^ 1 — > 
(see fig. [4]) , as opposed to the small negative value found 
for the corresponding solutions of Eq. ([2]). 

In principle, the results for F, D, and AF orientations 
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Figure 7: (Color online) General view of extrapolated phase 
diagram. Full line is extrapolation of solutions of Eq. (|2|), with 
the TM along D, using y ir = -1. Points: PRG, N = 14, with 
the TM along F (triangles) and D (hexagons). Inset: Low- 
temperature section of same data. Same axis labels and sym- 
bol captions as in main figure, except for additional dashed 
line which is extrapolation of solutions of Eq. (fj]), with the 
TM along F, using — — 1, 



must eventually extrapolate to the same location of the 
critical boundary as N — > oo. The way finite- A data 
vary as N increases is consistent with this, see Figs. Q] 
and [3] However, it is apparent that subdominant cor- 
rections to scaling have much larger amplitude for AF 
than for F or D. Attempting a theoretical understand- 
ing of why this is so usually becomes a highly nontrivial 
task, as it involves (i) unequivocally identifying the asso- 
ciated irrelevant exponents, and (ii) once this is done, 
analysing the (non-universal) amplitudes of the corre- 
sponding terms. Examples of this can be seen in Refs. [3 
and[l6|. Here, for practical reasons we chose to extrapo- 
late only the sets of (T, H) data for which the small sub- 
dominant corrections could be satisfactorily dealt with 
via Eq. |(4]) and its underlying assumptions. 

We have found that the minimum amount of discrep- 
ancy among all our results corresponds to the set of: (i) 
extrapolated data from the solutions of Eq. ([2]) , with the 
TM along D, using y ir = -1; (ii) PRG with the TM along 
D, and (iii) PRG with the TM along F. See Fig. The 
extrapolations from the solutions of Eq. ((2]) , with the TM 
along F, and using y ll: — —1, also agree very well with 
these, except for the low-temperature 'shoulder' where 
curves begin to depart more significantly from the hor- 
izontal line H = 2. This is illustrated in the inset of 
Fig. El 

In order to check on the universality of critical proper- 
ties, we examined the thermal exponent yr = Ijv along 
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Figure 8: (Color online) Thermal scaling exponent yr = 1/v 
along extrapolated critical line shown in FigO calculated via 
Eq. (©, with M = N - 2. TM along D. The discontinuity at 
T w 1.74 marks the change in assumed scaling direction (see 
text). 



the critical line. Considering two strips of widths M and 
N, finite-size scaling [l5[ gives: 



Vt = 1 



\u{k' n /k' m ) 
\n{M/N) 



(5) 



where 



V JV> 



"M 



are derivatives of the inverse corre- 



lation lengths, taken with respect to the appropriate 
temperature-like scaling field, evaluated on the critical 
curve. With the TM along D, and using M = N - 2, 
we swept the extrapolated line shown in Fig. ]7\ For sim- 
plicity, the temperature-like direction was taken as the 
temperature axis for low and middle fields H < 1.41 (cor- 
responding to 1.74 < T < 2.269 . . . ), and as the if-axis 
for the remainder of the critical line. The resulting so- 
lutions to Eq. ©, for TV = 10, 12, and 14, are shown 
in Fig. El Although the abrupt jump at T « 1.74 is an 
artifact, reflecting the above-mentioned (and somewhat 
arbitrary) change in the assumed scaling direction, one 
sees that on both sides of the discontinuity the estimates 



are rather close to the Ising value yx = 1, and systemati- 
cally approach it with increasing N. On the other hand, 
for T < 0.5 (where H already differs by less than 0.5% 
from the zero-temperature H c (0) = 2), crossover effects 
related to the energy level crossings at T — cause an 
extreme deterioration in our estimates. 

Near the H = extreme of the critical curve, we have 
fitted our extrapolated curves to a parabolic shape: 



T c {H) = T c (0)-aH 2 



(6) 



from which we get a — 0.217 ± 0.001, to be compared 
with a = 0.1767... and a = 0.3018..., each coming 
from a slightly different implementation of the interface 
method |7j. Note that the coefficient of a hypothetical 
linear term in Eq. ([6]) vanishes identically because the 
phase diagram is symmetric under field inversion, thus 
the scaling variable must depend on H 2 . For H — > 
this means that T C (H) = T c (0) — aH 2 plus higher-order 
terms [H]. 

Finally, we return to the exponentially vanishing de- 
viations from a horizontal line near T = 0, predicted in 
earlier work and mentioned in the Introduction. These 
are of the general form Q 



H = 2 + cT x exp(-d/T) 



(7) 



where c < for both implementations of the interface 
method 0], as well as for a free-fermion approxima- 
tion rjj] .whereas c > for the linear-chain approxima- 
tion [l|, lZJ (|c|, d, and x turn out to be of order unity in 
all cases). While we do not have enough accuracy at low 
temperatures to probe for this sort of effect, it seems safe 
to state that any stronger deviations from the horizontal, 
be they in the shape of a reentrance or the opposite, are 
ruled out by our results. 

Our extrapolated data for the location of the critical 
line, for both F and D, with y- 1T — —1 and —2, are avail- 
able as ASCII files 1191. 
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